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1 Introduction 

Although the main goal of the LHC machine is the search for and the study of the physics 
beyond the Standard Model, the expected huge rates of bottom quark production make es- 
pecially interesting the foreseen B physics programme for the LHC project. In fact a specific 
experiment (LHCb) will focus on B physics, while the two general-purpose experiments AT- 
LAS and CMS will dedicate special periods for data taking to this aim (see for example Ref. 

Among heavy flavour physics, heavy quarkonia production and decays have historically 
played a very important role in the development and test of Quantum Chromodynamics 
(QCD) as the best candidate to account for the strong interaction dynamics, and likely will 
continue keeping an outstanding position in this task. 

Moreover, over the last decade hadropro duct ion of heavy quarkonia has received a lot 
of attention from both theoretical and experimental viewpoints, to explain the discrepancy 
between the so-called colour-singlet model (CSM) and the experimental data, amounting 
to a factor of about 50 for direct J/ip hadroproduction at the Tevatron. In particular, 
the colour-octet mechanism (COM) || can be viewed as the (relativistic) generalization 
of the CSM and hence the most natural explanation for the unexpected surplus of heavy 
resonance hadroproduction. Nevertheless, when applied to other production processes like 
photoproduction at HERA, problems initially arose which cast doubts on the validity of 
the COM, although recent progress has been done allowing for a better understanding of 
the situation ||. Furthermore, results from Tevatron on charmonia polarization (one of the 
foremost predictions of the COM) seem to indicate even the failure of a naive application of 
the colour production mechanisms for charmonia jp. 

However, the COM can be viewed as deriving from a low energy effective theory, the Non 
Relativistic QCD (NRQCD) ||, so the question actually arising is whether NRQCD is the 
correct framework to deal with quarkonia production and decay. Perhaps the v expansion 
does not converge well for charmonium and subleading contributions cannot be neglected; 
perhaps the heavy quark spin symmetry is broken to a larger extent than expected. More 
work in this regard is required to clarify the situation. On the other hand possibly NRQCD 
is appropriate to describe bottomonia states and their production because of the larger mass 
of the bottom quark. Hence checking the COM in bottomonia hadroproduction is one of 
the challenges of strong interaction physics over the next years. Indeed there are alternative 
models in the literature (based on QCD) trying to explain the experimental facts (see for 
example || fjj). More astringent tests of heavy quarkonia production are thus required to 
enlight the situation, which can be qualified as rather confusing at present ||. 

In a series of previous works (||, || |ll], [l^, |13|) we have extensively analyzed char- 



monium hadroproduction in a Monte Carlo framework, using PYTHIA 5.7 |L4], |15| event 
generator with the colour-octet model implemented as a new routine in the generation code 
||10|| . Basically, such a production mechanism is based on the formation of an intermediate 
coloured state during the hard partonic interaction, evolving non-perturbatively into phys- 
ical heavy resonances in the final state with certain probabilities governed by NRQCD 0. 
In this work we extend our previous study of the T(1S) resonance [16|, [T^j to the whole 
T(nS) family below open bottom production, i.e. n = 1,2,3, using the CTEQ4L parton 
distribution function (PDF). A similar analysis can be found in Ref. [H| although limited 
to transverse momentum (py) values higher than 8 GeV. 

Although the discrepancies between the CSM and experimental cross sections on bot- 
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tomonia hadropro duct ion are smaller than those found for charmonia [19], still some extra 



contribution should be invoked to account for the surplus observed at the Fermilab Teva- 



tron. However, we find that, analogously to the charmonium case ||10|| , those matrix elements 
(MEs) determined from Tevatron data in other analyses |2(J have to be lowered once initial- 
state radiation of gluons is taken into account. This is because of the raise of an (effective) 
intrinsic momentum (kx) of the interacting partons enhancing the moderate and high-p^ 
tail of the differential cross section for heavy quarkonia production (for more details see 
Ref. ]|K|). This effect, as generated by the appropriate PYTHIA algorithm |TJ], ^TJ, is 
more pronounced - and likely more sound from a physical viewpoint - than a pure Gaussian 
smearing with a (required) large < kx > value. Besides PYTHIA, in smoothing the pro- 
duction cross section, endows us with the possibility of extending our analysis to the small 
Pt region of bottomonium production, keeping the assumption on the validity of the cross 
section factorization. 

The study of bottomonia production in hadron colliders should permit a stringent test of 
the colour-octet production mechanism, particularly regarding the predicted (mainly trans- 
verse) polarization of the resonance at high-pr HH, whereas other approaches, like the 
colour evaporation model, predict no net polarization; indeed, LHC experiments will cover 
a wider range of transverse momentum than at the Tevatron, allowing to explore the region 
p\ » 4mf , where m& denotes the bottom quark mass. 

In this paper we also present the prospects to probe the gluon density of protons via heavy 
quarkonia inclusive hadroproduction at high transverse momentum in the LHC. Our proposal 
should be viewed along with other related methods of constraining the gluon distribution in 
hadrons like di-jet, lepton pair and prompt photon production [E3|, [24], |25[. We must clearly 



state that it relies on the dominance of a particular production mechanism at high px (the 
COM) predicting a dominant contribution from gluon fragmentation. In spite of this and 
other assumptions (such as the validity of the factorization of the cross section), our feeling 
is that LHC collaborations should keep an open mind on all the possibilities offered by the 
machine, thereby exploring the feasibility of this proposal. For all these and other physical 
reasons, it is worth to estimate, as a first step, the foreseen production rate of bottomonium 
resonances at the LHC and this constitutes one of the goals of this work. 

We have based our analysis of bottomonia inclusive production on the results from Run 
IB of the CDF collaboration ][II| [26|. at the Fermilab Tevatron. This means significantly 
more statistics than the data sample from Run IA, employed in a former analysis P0fl . 
However, the different sources of prompt T(IS') production were not yet separated along the 
full accessible p^-range, in contrast to charmonium production. Hence we give in Section 2 
the numerical values for some relevant combinations of long-distance MEs, including direct 
and indirect T(nS) inclusive production, extracted from the fit to the CDF experimental 
points. (Prompt resonance production includes both direct and indirect channels, the latter 
referring to feeddown from higher T(nS) and Xbj{nP) states.) Nevertheless, we still are able 
to estimate some colour-octet MEs for direct T(IS') production from the measurements on 
different production sources at pt > 8 GeV [PTfl . The extrapolation to LHC is shown in 
Section 3 where we show the predicted differential and integrated cross section for all T(nS) 
resonances. In Section 4 we discuss heavy quarkonia inclusive hadroproduction as a probe of 
the gluon density in protons. Finally, in the appendices at the end of the paper, we gather 
those technical details and values of the parameters employed in the generation. 
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Figure 1: Different fits to the Tevatron data on bottomonia inclusive production in the rapidity 
interval \y\ < 0.4 using CTEQ4L PDF and m b = 4.88 GeV. First row panels: T(1S*); Second row 
panels: T(2S); Third row panels: Y(3S). Dot, dot-dash and solid lines correspond to the CSM, 
COM ( 3 S'^ only) and all contributions, respectively. The triangle mark indicates the pt lower cut- 
off used in the fit for each case: 2, 4 and 8 GeV. However, we plot the resulting curves extrapolating 
back over px > 1 GeV in all cases. 
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2 Fits to Tevatron data 



As briefly outlined in the Introduction, the theoretical differential cross sections on inclusive 
production of bottomonia would stand above Tevatron experimental points for relatively 
high pt if the set of long-distance parameters from |2(| were "blindly" employed in the 
PYTHIA generation running with initial-state radiation on. This is the analogous conclusion 
which one of us (M.A.S.L.) reached to in the equivalent analysis performed on charmonia 
hadroproduction |ID| . Indeed the smearing caused by multiple emission of gluons by the 
interacting partons is not limited to small pt values as could be initially thought, but its 
influence spreads over a larger region of transverse momenta. In fact we have checked, from 
a fit to the T(3S) differential cross section, that actually this effect amounts to a pretty 
large value for the effective < k F > of about 2 GeV, as we shall discuss more extensively in 
Section 2.2. 



Therefore we performed a new analysis of bottomonia CDF data [pq] , incorporating 
both direct and indirect production through the CSM (as a "fixed" contribution which, in 
fact, is dominant at low and even moderate pr) and the COM, adjusting the long distance 
parameters for different cut-offs from best \df=X 2 I^df fits to the experimental points, 
using the CTEQ4L PDF. If not explicitly stated the contrary, we are turning on initial- 
state radiation in addition to a Gaussian primordial kj- distribution (with a = 0.44 GeV, by 
default in PYTHIA) in all generations. 

2.1 Extraction of the colour-octet MEs 

In Figure 1 we show the theoretical curves obtained from our fits to CDF data. In general, 
nice fits, with \df values not too far from unity, were found, especially in the T(3S) case; 
instead, the Y(2S) came out to be the worst one. Let us stress that in the fitting procedure we 
excluded any possible negative contribution from the different channels at the cross section 



level, in contrast to [pL8| . Hence we had to dismiss any contribution from the So + Pj 
channels or, in other words, we set the M5 long-distance parameter (as defined in appendix 
A) equal to zero, since any positive contribution from this channel would lead to a worse 
X DF value in all cases. 

Table 1: Values of < 0j (n5) ( 3 Si) > \ tat ; n = 1,2,3 (in units of 1(T 3 GeV 3 ) from the best fits to 
CDF data at the Tevatron on prompt T(nS) inclusive production for different px lower cuts. We 
also provide the Xdf value for each case. The CTEQ4L PDF was used with initial-state and AP 
evolution activated in PYTHIA. 



Pt cut-off: 


2 GeV 


Xdf 


4 GeV 


Xdf 


8 GeV 


Xdf 


IS 


77±17 


1.74 


87±16 


1.53 


106±13 


1.00 


2S 


40±29 


2.87 


73±18 


1.58 


103±27 


1.87 


3S 


99±11 


1.00 


91±15 


1.00 


68±11 


1.00 



In Table 1 we show the values of < 8 ( 3 £i) > \ tot {n = 1,2,3), as defined in (A. 7), 
for different p? lower cut-offs, in correspondence with the plots of Figure 1. All values are 
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roughly of order 1CF 1 GeV 3 and agree, within the errors, with the results obtained for pt > 8 
GeV by the authors of Ref. |L8| using the CTEQ5L parton distribution function. 

Nevertheless, let us stress that our numerical estimates for the colour-octet MEs have to 
be viewed with some caution because of the theoretical and "technical" (due to the Monte 
Carlo assumptions) uncertainties. For example, our algorithm for AP evolution (see appendix 
B) should be regarded as a way of reasonably steepening the high-p^ tail of the (leading- 
order) differential cross section, which otherwise would fall off too slowly as a function of 
Pt- 



2.1.1 Separated production sources for px > 8 GeV 

Current statistics does not permit to subtract indirect production sources to obtain the direct 
T(IS') production cross section along the full accessible p^-range. Nevertheless, feeddown 
from higher states (T(nS'), Xbj( n P)) was experimentally separated out for pt > 8 GeV 
fi~9l , 27|. We used this information to check our analysis a posteriori (rather than using it 



as a constraint in the generation) and to draw some important physical conclusions. To this 
end the relative fractions of the contributing channels for pt > 8 GeV are reproduced in 
Table 2 from Ref. 



|19|, g7| . On the other hand, we show in Table 3 (which updates our results 
I~6|j using CTEQ2L) the fractions found in this work corresponding to the 
different generated channels for p? > 8 GeV, following the notation introduced in appendix 
A. 



presented in Ref. 



Table 2: Relative fractions (in %) of the different contributions to T(15*) production from CDF 



data at pt > 8 GeV [27|. Statistical and systematic errors have been summed quadratically. 



contribution 


Tevatron results 


direct T(1S) 


50.9±12.2 


T(25)+T(3S) 


11.5±9.1 




27.1±8.2 


Xb(2P) 


10.5±4.6 



Table 3: Relative fractions (in %) of the different contributions to Y(IS') production at the Tevatron 
for px > 8 GeV from our generation. Possible contributions from Xbj(3P) states were not generated. 



contribution 


our generation 


T(IS') 3g (8) 


36.8 


T(1S)I CS M 


19.5 


T(2S)+T(3S)\ CSM 


3.9 


Xb(lP)\ CSM 


24.1 


Xb{2P) csm 


15.7 



By comparison between Tables 2 and 3 we can conclude that the T(IS') indirect pro- 
duction from Xfoj's decays is almost completely accounted for by the CSM according to the 
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assumptions and values of the parameters presented in appendix A. Indeed, experimentally 
37.6±9.4% of Y(IS') production is due to Xbj(lP) an d Xw(2-P) decays while from our 



generation we find a similar global value, namely 39.8%, coming exclusively from colour- 
singlet production! Moreover, assuming that a 7.6% from the 36.8% fraction (corresponding 

to the colour-octet 3 S[ 8 ^ contribution as expressed in Eq. (A. 7)) can be attributed to the 
T(2S) + T(3S) channel in addition to the colour- singlet contribution (3.9%), we obviously 
get the fraction 11.5% for the latter, bringing our theoretical result into agreement with the 
experimental value. This single assignment implies to reproduce quite well the experimental 
fraction (« 51%) of direct T(IS') production by adding the remaining 3 s[ 8S) contribution to 
the T(1S) CSM channels (~ 49%) in our generation. 

Of course all the above counting was based on mean values from Table 2 and subject 
to rather large uncertainties. Nevertheless, apart from the consistency of our generation 



w.r.t. experimental results under minimal assumptions, we can conclude again as in |16 
that there is almost no need for T(IS') indirect production from feeddown of XbJ states 
produced through the colour-octet mechanism. In other words, the relative contribution 
from P-wave states to < Oj^ 1S \ 3 Si) > \ tot in Eq. (A. 7) should be quite smaller than naively 
expected from NRQCD scaling rules compared to the charmonium sector, in agreement with 
some remarks made in p8| and recent results found in |18|]. The underlying reason for this 



discrepancy w.r.t. other analyses |2(| can be traced back to the dominant colour-singlet 



contribution to the cross section at values as much large as ~ 18 GeV (see Figure 1) 
caused by the effective smearing - already applied to charmonium hadropro duct ion by 
one of us ||T0 . 



On the other hand the corresponding velocity scaling rule in the bottomonium sector is 
roughly verified as we shall see. Defining the ratios of matrix elements: 

Rv(n) = ° J l} Lt , (1) 
< O^SJ > \ tot 

its values, shown in Table 4, are in accordance with the expected order-of-magnitude v 4 ~ 0.01, 
where v is the relative velocity of the bottom quark inside bottomonium. Nevertheless we 
realize an increase of R v (n) for higher n values. Assuming that the < Oj^ nS \ 3 Si) > \ tot 
matrix element could be interpreted as a (weighted) colour-octet wave function squared (in 

the same way as < Oj^^^Si) > \ tot w.r.t. the colour- singlet state, see appendix A) the 
ratio R v (n) of both squared wave functions in the origin comes out as not independent of 
the resonance state under consideration. 



Table 4: Values (in units of GeV 3 ) of different colour-singlet and colour-octet combinations of MEs 
according to Eqs. (A. 4) and (A. 7) and the ratios R v (n); n = 1,2,3. The best Xdf values from 
Table 1 are displayed. 



Resonance 


<or»\*s 1 )>\ tot 


<O^^S 1 )>\ tot 


R v (n) 


T(1S) 


11.1 


0.106 


0.0095 


T(2S) 


5.01 


0.073 


0.0145 


T(3S) 


3.54 


0.099 


0.028 
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Thus we conclude that this particular NRQCD velocity scaling rule, although valid as an 
order-of-magnitude estimate, retains a weak dependence on the principal quantum number 
n, not completely cancelling in the ratio (1). 



T(3S) CTEQ4L (AP) kt= 1 .8 GeV 




Figure 2: Fit to the Tevatron data on T(35 l ) hadroproduction using a gaussian smearing function 
with a = 2 GeV, i.e. < k T >= 1.8 GeV. 



2.2 Gaussian < kx > smearing 

The smearing effect on the differential cross section caused by initial-state radiation of gluons 
can be roughly simulated by means of a gaussian intrinsic kx distribution of the interact- 
ing partons inside hadrons, to be convoluted with the corresponding hard interaction cross 
sections: 

D(k T ) = -Lexp(-§) (2) 

with 



< h: > = ^ o (3) 



The width of the gaussian can be viewed as an adjustable parameter |2J| . In fact PYTHIA 
incorporates as an option a gaussian primordial kx smearing, whose width can be set by the 
user. We used this possibility to make a "new" fit of Tevatron data for the T(3S) resonance, 
employing the same matrix elements as shown in Table 1 but with initial-state radiation 
off. Then the gaussian k^ smearing has to simulate the (this time missing) initial-state 
radiation. In Figure 2 we show the resulting histogram, corresponding to a value a = 2 
GeV, i.e. < k T >= 1.8 GeV. 



3 T(nS) Production at the LHC 

We already mentioned in the Introduction that bottomonium hadroproduction is especially 
interesting to check the validity of the colour-octet model as often emphasized in the litera- 
ture |30|. 31 1 . This becomes particularly clear at the LHC since experimental data will spread 
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Figure 3: First row panels, left side : Predicted prompt T(15) differential cross section (multiplied 
by the muonic branching fraction) at the LHC using the CTEQ4L PDF and m& = 4.88 GeV. A 
rapidity cut \y\ < 2.5 was required for bottomonium. Dot-dashed line: contribution; solid 

line: all contributions. Right side : Integrated cross section. Second and third row panels: The 
same as in first row for T(2S) and T(3S) respectively. 
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over a wider p^-range than at the Tevatron, allowing an overall study from low to very high 
Pt values. Therefore the expected transition of the different production mechanisms along 
the pt region could be scrutinized in detail: from gluon gluon fusion at low pt to the foreseen 
asymptotically dominant gluon fragmentation into bottomonium states. 

Keeping this interest in mind, we used our code implemented in PYTHIA to generate 
prompt T(nS) resonances in proton-proton collisions at a center-of-mass energy of 14 TeV 
employing the best Xdf colour-octet MEs shown in Table 1. In figure 3 the theoretical 
curves for the T(nS) (n = 1, 2, 3) differential and integrated cross sections are exhibited as a 
function of px, including both direct production and feeddown from higher resonance states 
(except for the T(3S)). 

In Figures 4 we show our prediction for direct T(nS) production. This is especially 
interesting if LHC detectors would be able to discriminate among those different sources of 
resonance production. (See the end of Section 4 and footnote #2.) 

To this end we generated T(IS') events through both the CSM and COM making use of 
the following parameters 

• < Oj {1S \ 3 S 1 ) > \ direct = 9.28 GeV 3 (from @) 

. <0^ ls \*S 1 )>\ dtrect = 0.084 GeV 3 

The first value corresponds to the CSM matrix element for direct production while the 
< Oj < - 15 ^( 3 S'i) > \ dir value was obtained after removing the T(2S) + T(35) contribution 
according to the discussion made in Section 2.1.1, i.e. under the assumption that a fraction 
7.6% from the 36.8% in table 3 should be assigned to indirect production. Finally let us 
mention that we neglected any contribution from the 1 Sq 8 * > + 3 Pj 8 ^ channels, in accordance 
with our analysis on Tevatron results of Section 2. 




Figure 4: The same as in Figure 3 for direct T(15") production at the LHC. 
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4 Heavy quarkonia inclusive production as a probe of the gluon 
density in protons 

One of the goals of the LHC project is to perform precise tests of the Standard Model of 
strong, weak and electromagnetic interactions and the fundamental constituents of matter. 
In fact the LHC machine can be viewed as a gluon-gluon collider to some extent. On the 
other hand, many signatures (and their backgrounds) of physics, both within and beyond 
the Standard Model, involve gluons in the initial state |]J . Therefore an accurate knowledge 
of the gluon density in protons acquires a special relevance for all these reasons. 

So far, the most precise determinations of the gluon momentum distributions in the pro- 
ton come from the analysis of the scaling violations of the structure function F2. However, 
this represents an indirect method since it is the sea distribution which is actually measured 
and the gluon density is obtained by means of the QCD evolution equations. On the other 
hand, hadron-hadron scattering processes with direct photon production or jets in the fi- 
nal state will probably be extremely adequate to probe "directly" the gluon distribution in 
hadrons. In this Section, we shall examine the possibility of using heavy quarkonia inclusive 
production in proton-proton collisions at the LHC, in a complementary way to those stud- 
ies. However, there are still many uncertainties and pending questions regarding quarkonia 
hadroproduction despite the existence of NRQCD ||, an effective theory coming from first 
principles, which should provide an adequate framework for this kind of processes involving 
both perturbative and non-perturbative aspects of the strong interaction dynamics. Likely, 
forthcoming experimental data - and their respective analyses - from Tevatron and other 
machines like HERA, should clarify the situation before LHC becomes operative. 

In the following we shall focus on very high production of bottomonia states. Therefore 
the main production mechanism according to the COM should be through the partonic 
subprocess: 

9 9^ 9* 9 (4) 
followed by the gluon fragmentation into a T(nS) state: 

g*^T(nS)X ; (n = l,2,3) (5) 

produced through a colour-octet mechanism. On the other hand, the bottom mass is large 
enough to justify the colour-octet model applied to quarkonium hadroproduction. Whether a 
similar approach could be applied to charmonium resonances has to be checked, for example 
analyzing the transverse polarization of the resonance. 

Ideally, the final state gluon (g) in Eq. (4) will give rise to a recoiling jet (g^jet), 
sharing, in principle, the same transverse momentum as the heavy resonance (in absence of 
higher order corrections; see however appendix C). Hence events would topologically consist 
of an almost isolated muon pair from the decay of the heavy resonance and a recoiling jet. 
Indeed one should expect a pair almost isolated because the energy difference between 

the masses of the intermediate coloured and final states is assumed to be rather small (of 
the order of m^v 2 ~ 500 MeV) then allowing the emission of eventually a few light hadrons 
via soft gluon radiation at the final hadronization stage. 

Bottomonia production coming from fragmenting gluons in QCD jets (an alternative 
production mechanism, see P2[ ) should not exactly display the same signature as the hard 
o? s processes. Indeed, the muon pair would be embedded in one of the two jets - not so much 
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isolated as in the process (4-5) due to the production cascade - and its momentum should 
not balance the momentum of the other event jet to the same extent. In sum, the signature 
of an almost isolated muon pair recoiling against a jet with an approximate momentum 
balance in the transverse plane, should provide a suitable tag for the production mechanism 
represented in Eqs. (4-5). 

We shall assume a tight kinematic cut in this our first approach: Both the rapidity 
of the heavy resonance and the rapidity associated to the recoiling jet should be around 
zero. However, in order to increase the foreseen statistics, one could dispense with this 
constraint by only requiring (within the experimental and theoretical uncertainties) back- 
to-back production. We shall come back to this issue at the end of this Section. 



Developing the idea 

In the absence of any intrinsic kx effect, we can write the triple differential cross section for 
the inclusive production process pp-^TX as 

j ° = 2p T J2XaXbfa/p(x a )f b/p (x b )^^ (6) 

dy r dy jet dp T ^ dt 

where f a / P (x a ) denotes the parton-a density in the proton, and 



da ab da 1 
dt dt V ' 16tts 2 



^(^ Tc ) = T7-^El^-^)| 2 (7) 



stands for the partonic differential cross section (the barred summation denotes an average 
over initial and final spins and colours) consisting of a short distance (and calculable) part 
and a long distance part which can be identified as a colour-octet matrix element according 
to NRQCD. This factorization of the cross section was established on solid grounds in Ref. 
H within the NRQCD framework. 

As above-mentioned we shall require both rapidities (of the T and the recoiling jet) to 
be less than a common small value y Q : \yy\ < y , \yj e t\ < Vo- We could set y = 0.25 for 
example, as discussed in appendix C.) Then x a ~ x b = x, and 



.,2 



S 



x a x b X 

s 

At very high pt (i.e. p\ » Ami) we can identify s m A.p\. (Hereafter we consider pt > 20 
GeV.) Therefore measuring the transverse momentum of the resonance should lead to the the 
knowledge of the momentum fraction x of the interacting partons, with a typical uncertainty 
(see appendix C) 

Ax , , 

— = yo 9 

x 

In particular the dominant partonic subprocess should be the gluon-gluon interaction. 
Thus the gluon density G(x, /i 2 ) in the proton will mainly be involved and we can write as 
a first approximation 

. f\ =2 Pr x 2 (W) 2 % (10) 
dy r dy je tdp T dt 

where we can choose, for example /i 2 = s 
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The proposal [] 

We propose to study the ratios: 

x\ G(x 2 ,fij) 2 _ / dcr gg / dti \ ^ /pri \ / d 3 a/ dy r dy jet dp T 2 \ ^ 
x 2 G(a?i, /if) 2 \da g g/di 2 J \PtzJ \d 3 a /dy r dy jet dp T iJ 

for a set of £i,:E2 pairs and different gluon distributions. The number of pairs is basically 
limited by Ax, i.e. y , so this constraint cannot released too much (see appendix C). 

Therefore the keypoint is to consider the l.h.s. of the above equality (Eq. (11)) as an input 
corresponding to different sets of the gluon distribution for the proton, whose x dependence 
is hence assumed to be "known", and in fact would be tested. On the other hand the r.h.s. 
corresponds to an input from experimental data and some theoretical factors likely under 
control. 

Let us remark that the x and fi 2 values are not independent in this proposal; indeed for 
each value of x, /i 2 is fixed by s = x 2 s. However, notice that the scale can actually be varied 
by choosing a different assignment for /z 2 , e.g. fi 2 = s/4. 

Next we shall write expression (11) as 



2\2 



x 2 G(x 2 ,f4) 



x\ G(x 1 ,fi 2 1 ) 2 

where 

d&ggjdt 



Rtheo x Rexp (12) 



Rtheo\JPtUPt2, fh, Hz) = /corX- — j- (13) 

aaggjat 2 



and in the high p? limit, 

Rtheo ~^ f, 



««(^D PT2 



a 3 sit4) Pti 

explicitly showing that a s (fi 2 ) is entangled in the gluon density determination. On the other 
hand, note that the dependence on the NRQCD matrix elements does cancel in Rtheo, but 
there is a dependence on the scales /x 2 and which should match the same dependence in 
the left hand side. We have incorporated some possible corrections through the f cor factor 
- which could be calculated either analytically or by Monte Carlo methods - taking into 
account higher-order effects such as intrinsic of the interacting gluons, AP evolution of 
the fragmenting gluons, etc. 

On the other hand the experimental input reads as the ratio 

o / \ (Ptx\ v fd 3 a/dyrdy je tdp T2 \ , 

Rexp{PTl,PT2, yo)=[ X / (14) 

\PT2 J \d 6 a/ dy T dy jet dp T i J 
which can be obtained directly from experimental data. 

Tresented at the UK Phenomenology Workshop on Heavy Flavour and CP violation, Durham, September 2000 
[p3| and at the B physics working group meeting of the ATLAS collaboration held at CERN in October, 2000 
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Introducing the gluon quark contribution 

Although expectedly dominant, the gluon gluon partonic subprocess is not the only 
contribution to the cross section yielding a fragmenting gluon into T(nS) at high p?. Also 
gluon quark scattering gq-^g*q followed by g*— >T(nS)X, can give a sizeable contribution 
(about 20% at pr > 20 GeV, see table 6 in appendix A). Consequently, the expression (12) 
for the ratio of gluon densities has to be modified to include the quark distribution q(x,fj, 2 ) 
in protons: 

x 2 G{x 2 ,fi 2 ) {x 2 G{x 2 , fi 2 2 ) + k-x 2 q{x 2 , f4)) R x R , 15 x 
XiG(xi,/xf) (x 1 G(x 1 ,/i 2 ) + k-x 1 q(x 1 ,nl)) the ° exp 

where k is a factor taking into account the ratio of the gq and gg cross sections, both 
calculated at the same values of the Mandelstam variables s and t of the hard interaction, 
i.e. 

k = (16) 

d(Jgg/dt 

becoming independent of x (and /i 2 ) at zero rapidity and large pr] then k ~ 0.2. 
Alternatively, one can write the density ratio as 

x\ G{x 2 , n 2 2 ) 2 {l + k-X{x 2 , nl)) 

xlGixuWil + k-Xixurt)) { } 

where 

G(x, n z ) 

By Taylor expanding the above ratio, the leading term is 

x\ G(x 2 ,f4) 2 
x\ G(x 1 , /if) 2 

where 



;i + r) 



r = kx 



\(x 2 ,i4) - A(xi,/i 2 ) 



should be a quite small quantity. (We have checked with CTEQ4L that typically r w 0.03 
for values between x 1 = 3-10 -3 and x 2 = 1.5-10 -2 .) 
We can rewrite Eq. (12) as 



x\ G(x 2 ,n 2 2 ) 2 



x 2 G( Xl , fij) 2 " 1 + r) ~ Rthe ° X Rexp (18) 

Again the l.h.s. is an input from the PDF to be tested, while the r.h.s. comes from experi- 
mental data and some theoretical calculations without requiring the NRQCD MEs values. 

From an experimental point of view it may happen that the discrimination among the 
different T(nS) states via mass reconstruction could become a difficult task, especially at 
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very high p?, because of the uncertainty on the measurement of the muons momenta (lj f\ 
Nevertheless, since we are proposing to study ratios of cross sections, we can consider the 
overall T(nS) inclusive production, without separating the different bottomonia sources - 
all the weighted matrix element cancelling in the quotient if we neglect the mass differences 
between the different states. (Notice that at high px there is almost no contribution from 
the CSM.) In Figure 5 we show the combined production rate at p? > 20 GeV and \y\ < 0.25 
for the upper and lower values of the colour-octet matrix elements shown in Table 1. 




Figure 5: Predicted T(1S) + T(25) + T(3<S) weighted contributions to bottomonia inclusive produc- 
tion at the LHC corresponding to the upper and lower MEs from Table 1, in the rapidity interval 
\y\ < 0.25 and px > 20 GeV. Left panel: differential cross section; right panel: integrated cross 
section. 

Assuming an integrated luminosity of 10 fb _1 , corresponding to one year running (10 7 s) 
of LHC at "low" luminosity (10 33 cm~ 2 s _1 ) we can easily get the expected number of events 
from figures 5, just by multiplying the ordinate by a factor 10 7 . Thus we can see that the 
foreseen number of events (aside efficiency reduction) at px > 20 GeV is about 10 5 , whereas 
at pt > 40 GeV is about 10 4 . By extrapolation we get a meagre expected number of ~ 10 2 
events at pt > 100 GeV. This makes unlikely any measurement for transverse momentum 
larger than 100 GeV, under the tight rapidity cut of 0.25 on the resonance which we are 
imposing. 

In view of the foreseen rates of bottomonia production at the LHC we propose testing the 
shape of the gluon density in protons for x values ranging in the interval: 3 x 10~ 3 to 1 .5 x 10 -2 , 
using x = Js/ s from p T = 20 GeV up to px = 100 GeV, under the rapidity constraint 
y < 0.25. 

Nevertheless, by removing the condition \y\ < 0.25 statistics could considerably be en- 
larged. Since our proposal essentially relies on the determination of the Feynman x of the 
interacting partons by measuring the pt of the final products of the reaction, there is still the 
possibility of requiring a back-to-back topology but sweeping the whole accessible rapidity 
region \y\ < 2.5, instead of limiting ourselves to the central rapidity values. This goal can 
be achieved by selecting events with the muon pair and the recoiling jet sharing common 

2 We thank M. Smizanska for bringing this point to our attention. 
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values of pr and absolute rapidities, within the uncertainties. In other words, events could 
be accepted with both T(nS) and recoiling jet rapidities satisfying \ \y~c\ — \Vjet\\ < 0.25; in 
such a way statistics should increase by a factor ~ 10, possibly extending the allowed region 
of px up to higher values than 100 GeV, and hence reaching larger values of the momentum 
fraction x. 

As a final remark, if the colour-octet model is confirmed and the corresponding MEs 
accurately and consistently extracted from other experiments like Tevatron or HERA - or 
theoretically computed - one can consider then the possibility of unfolding the gluon density 
from the measured cross section as proposed, for instance, in Ref. |35|] by means of D* 
meson production at HERA. In such a case, our proposal would extend beyond the study of 
ratios, allowing the extraction of gluon and quark densities directly from heavy quarkonia 
production mechanisms. 

5 Summary 

In this paper we have analyzed CDF measurements on T(nS) {n = 1, 2, 3) inclusive hadropro- 
duction cross sections at the Tevatron in a Monte Carlo framework, extracting some relevant 
colour-octet NRQCD matrix elements. Higher-order QCD effects such as initial-state radi- 
ation and AP evolution of gluons were included in our analysis. The fact that we were able 
to take into account the effects of soft-gluon emission in the initial-state (according to the 
PYTHIA machinery) yielding a smooth pt behaviour at low p?, allowed us to extend the 
study of the differential cross section down to pt = 2 GeV. However, in order to study the 
sensitivity to different cut-offs, fits were performed using experimental points above several 
Pt cut-offs (see Table 1). 

On the other hand, since the different sources of T(nS) production were not experi- 
mentally separated along the full accessible pr-range we have included all of them in the 
generation and later fits. Only for pt > 8 GeV, feeddown from XbJ states was experimentally 
separated out from direct T(IS') production. We used these results as a consistency check 
of our analysis and to draw some conclusions summarized below. 

The numerical value of the < Og^ ( 3 5i) > \ tot matrix element should be ascribed almost 
totally to T(nS) states. This finding may be surprising when confronted with other analyses 
f20| , p8| , where the contribution to the T(IS') yield through the colour-octet Xbj channels 
was thought as dominant p8l [HL |36|. On the contrary, we concluded from Tables 2 and 3 
that the colour-singlet production can account by itself for the feeddown of T(1S) from Xbj 
states. (Notice however that experimental uncertainties still leave some room for a possible 
COM contribution but to a much lesser extent than previously foreseen pBfl.) On the 
other hand the different production channels are consistent (or can be made consistent) with 
the experimental relative fractions shown in Table 2, after some reasonable assumptions. 

We have extended our study to LHC collider experiments (y/s = 14 TeV center-of-mass 
energy). In Figure 3 we presented our predictions for prompt production rates (i.e. including 
direct and indirect production) while in Figure 4 we showed our prediction for direct T(1S) 
production alone. 

We conclude that the foreseen yield of T(nS')'s at LHC energy will be large enough, 
even at high-p^, to perform a detailed analysis of the colour-octet production mechanism 
and should be included in the B-physics programme of the LHC experiments, probably 
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deserving (together with charmonia) a dedicated data-taking trigger. 

On the other hand, we have presented in some detail (but without considering detector 
effects) the prospects to use heavy quarkonia inclusive production at the LHC (or perhaps 
at the Tevatron too, using data from the high luminosity Run II) to probe the gluon density 
in protons. Let us remark that we are renouncing to use any absolute normalization because 
this would imply a precise knowledge of the colour-octet NRQCD matrix elements governing 
the transition of the fragmenting gluon into a particular heavy quarkonium state, hence 
introducing an important uncertainty since they are not accurately known so far. If those 
NRQCD elements were finally accurately and consistently determined (either theoretically 
or experimentally) one could then consider the possibility of unfolding the parton densities 
from LHC experimental data on heavy resonance inclusive production, as described in our 
proposal. Thus we conclude that: a) Inclusive hadroproduction of heavy resonances at high 
Pt in general-purpose LHC experiments could be a complementary method of constraining 
the gluon density in protons, along with other related methods: di-jet, lepton pair and 
prompt photon production; b) An experimental advantage of this method lies in the fact 
that the "flight direction" of a high-pr fragmenting gluon into a T(nS) resonance can be 
inferred from the muonic pair coming from its decay, thereby providing a clean constraint 
in the search for the associated recoiling jet, in addition to the self-triggering signature of 
events. 

We finally want to stress the importance of keeping an open mind on the different pos- 
sibilities offered by the LHC, likely deserving a previous work to prepare in advance jointly 
experimental strategies and theoretical calculations. 
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Appendices 



Some technical details for the generation with PYTHIA 



Basically we reproduce here the information given in [16fl about the values of 
the parameters and options set in running PYTHIA, with some more details. 

Originally the event generator PYTHIA 5.7 produces direct J/ip and higher Xcj 
resonances via the CSM only fl4j| . It is not a difficult task to extend this gen- 



eration to the bottomonium family by redefining the resonance mass and wave 
function parameter accordingly. In our analysis we have besides implemented a 
code in the event generator to account for the colour-octet production mecha- 
nism via the following a 3 partonic processes: 

(A.l) 
(A.2) 
(A.3) 

where (QQ) [ 2S+1 Xj] stands for a certain heavy quarkonium state denoted by 
its spectroscopic notation (see Refs. ^(| [T(J for more details). In particular 
we have considered the 3 S[ 8 \ 1 5 , q 8 ' ) and 3 Pj contributions as leading-order 
intermediate coloured states. In addition we generated T(nS) [n = 1,2,3) and 
Xbj( n P) [n = 1,2) resonances decaying into T(IS'), according to the CSM as 
mentioned above. 



9 


+ 


9 ~ 


- (QQ)? S+l Xj] 


+ 9 


9 


+ 


Q ~ 


-> (QQ)[ 2S+l Xj] 


+ q 


Q 


+ 


1 - 


-> (QQ)[ 2S+1 Xj] 


+ 9 



Table 5: 3 S[ contributions to the T(15') cross section at the Tevatron for > 8 GeV 



Contribution 


% 


99 


69 


qg 


30 



qq 1 



Table 6: 3 S[ contributions to the T(1S) cross section at the LHC for > 8 GeV 



Contribution 


% 


99 


80 


qg 


20 


qq 


~ 



20 




», (GeV) 



Figure 6: Gluon-gluon versus quark-gluon 3 S[ contributions from our T(IS') generation at the 
Tevatron for pr > 8 GeV. The latter becomes more and more important at larger as could be 
expected since higher Feynman x of partons in the proton are involved and the gq contribution 
becomes increasingly more significant w.r.t. the gg one. 



A px lower cut-off was set equal to 1 GeV (by default in PYTHIA) throughout 
the generation since some of the contributing channels are singular at van- 



ishing transverse momentum [fCq| . Furthermore, all fits of Tevatron data were 
performed using p T values above 2 GeV. 

We find from our simulation (see Tables 5 and 6) that gluon-gluon scattering 
actually stands for the dominant process at high pt as expected, gluon-quark 
scattering contributes appreciably however (~ 20 — 30% of the colour-octet pro- 
duction cross section) whereas the quark- ant iquark scattering represents a quite 
small fraction ( ~ 1% at the Tevatron). In Figure 6 we plot the gluon-gluon 
and quark-gluon 3 S[ 8 ^ contributions as a function of the transverse momentum 
of the resonance obtained from our generation for the Tevatron. This kind of 
information could be particularly interesting for our discussion on the probe of 
the gluon density in protons developed in Section 4. 



Set of "fixed" and free parameters used in the generation 

Below we list the main parameters, including masses and branching fractions, 
used in our generation with PYTHIA 5.7. We employed the CTEQ4L parton 
distribution function (PDF) in all our analysis. 



Masses and branching fractions: 

• m b = 4.88 GeV 

• ^resonance — ^ m b m the theoretical calculation of the cross sections for 
the short-distance processes [p0|. In the phase space factors for the event 
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generation we set, however, real masses of the T(nS) and weighted mean 



values for the Xb(nP) resonances [16] . 
BR[r(lS)^fM + fM-] = 2.48 % (|3) 
BR[r(2S)^n+fi-] = 1.31 % 



BR[?(3S)^fj, + n-] = 1.81 % (0) 



Colour- singlet parameters (from [ZB|]): 

. <0 1 T(1S) ( 3 S 1 )>L = H-lGeV 3 
• < OJ (2B) ( 3 Si) > L, = 5.01 GeV 3 



I tot 



<0 1 T(3S) ( 3 S 1 ) > 



to* 



3.54 GeV 3 , defined as 



Lt = E < 0? (mS) ( 8 £i) > MT(mS)->T(rtf)X] 

m>n 



<Of 1(lp) ( 3 P 1 ) >=6.09 GeV 5 
<Of 1(2P) ( 3 P 1 ) >= 7.10 GeV 5 



(A.4) 



The radial wave functions at the origin (and their derivatives) used in the 
generation can be related to the above matrix elements as 

<Oj inS) ( i S 1 )>= ^-\R n (0)\ 2 (A.5) 



< o\ 



XbJ(nP) rtp\ > 



9 

2^ 



(2J+1)|<(0)| S 



(A.6) 



whose numerical values were obtained from a Buchmuller-Tye potential model 



tabulated in Ref. 38 



Colour-octet long-distance parameters to be extracted from the fit: 



< Oj (nS \ 3 S 1 ) > \ tot , defined as 



< Oj {nS \ 3 S 1 ) > 



T(mS)/3 



tot 



5i) > Br[T(mS)^T(nS)X] 



E <o 

m>n 

E E < Or AmP) ( 3 S 1 ) > Br[ Xb j(mP)^T(nS)X] 

m>n J=0 

(A.7) 



<0 8 T(n5) (%)> 



tot' 

v -So) > I tot 



defined as 

3 



< O^ {mS) CS ) > Br[T(mS)^T(nS)X] 

(a.s; 
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< Oj (nS \ 3 P ) > \ toV defined as 
< O 8 T( " 5) ( 3 P ) > \ tot = E < Oj {mS \ 3 P ) > Br[T(mS)^T(nS)X] 

m>n 

(A.9) 



On the other hand, the differences in shape between the 1 S < if l and 3 -Pj 8 ' ) contri- 
butions were not sufficiently great to justify independent generations for them. 
In fact, temporarily setting < Oj^ 1S \ 3 Po) >= ml < oJ^^Sq) > and defining 
the ratio 

■^2 da [3d (8)1 
dp T I D J 

it is found r ~ 5 as a mean value over the [0, 20] GeV p^-range. Actually the 
above ratio is not steady as a function of the T(IS') transverse momentum. 
Therefore in the generation we split the pt region into two domains: for p? < 6 
GeV we set r = 6 whereas for p? > 6 GeV we set r = 4. 

In summary, only the 1 S , q 8 ' ) channel was generated but rescaled by the factor r to 
incorporate the 3 -Pj 8 ' ) contribution as we did in [TO] for charmonium hadropro- 



duction. Consequently, in analogy to |2U| we shall consider only the combination 
of the colour-octet matrix elements: 

M5 = 5 x (<^" S, (%)>U + <QWo)>LA (A . n) 
V 5 J 
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CTEQ4L versus CTEQ2L 

In our previous work [|l6j on Y(IS') hadroproduction we employed the (now 
outdated) CTEQ2L parton distribution function, which however is still often 
used in current applications of PYTHIA at LHC collaborations. Throughout 
this paper we have used CTEQ4L but a comparison with the previous analysis 
is in order. In Figure 7 we present the two fits to the same Tevatron data [^6| . 
Notice that the 1 5'o + 3 Pj contribution has been disregarded in the latter case. 



Table 7: Colour-octet matrix elements (in units of 10 -3 GeV 3 ) from the best fits to CDF data at 
the Tevatron on prompt Y(IS') inclusive production, using either CTEQ2L [15] or CTEQ4L parton 
distribution functions respectively. 



ME: 


<0^ s \3 Sl )>\ tot 


M T ^-5x(< O " i T P0)> + <O « U " (l50)> ) 


CTEQ2L 


139±18 


6±5 


CTEQ4L 


77±17 


~ 



The CTEQ4L PDF incorporates a BFKL style rise at small x, rather than a 
flat shape as in CTEQ2L. Therefore it is not surprising that we find smaller 
values for the colour-octet matrix elements in the former case, as can be seen 
in Table 7. 
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B 

Altarelli-Parisi evolution 

According to the colour-octet model, gluon fragmentation becomes the domi- 
nant source of heavy quarkonium direct production at high transverse momen- 
tum. On the other hand, Altarelli-Parisi (AP) evolution of the splitting gluon 
into (QQ) produces a depletion of its momentum and has to be properly taken 
into account. If not so, the resulting long-distance parameter for the 3 S[ 8 ^ chan- 



nel would be underestimated from the fit fT2|j . 

The key idea is that the AP evolution of the fragmenting gluon is performed 
from the evolution of the gluonic partner of quarkonium in the final-state of the 
production channel 

9 + 9 - 9 m (-+(QWsF\) + 9 (B.l) 



Let us remark that, in fact, g* is not generated in our code |T(|. Final hadroniza 



tion into a (QQ) bound state is taken into account by means of the colour-octet 



matrix elements multiplying the respective short-distance cross sections [^0 , K| . 
Nevertheless, it is reasonable to assume that, on the average, the virtual g* 
should evolve at high px similarly to the other final-state gluon - which actu- 
ally is evolved by the PYTHIA machinery. We used this fact to simulate the 
(expected) evolution of the (ungenerated) g* whose momentum was assumed to 
coincide with that of the resonance (neglecting the effect of emission/absorption 
of soft gluons by the intermediate coloured state bleeding off colour |TT|). 

Therefore, event by event we get a correcting factor to be applied to the trans- 
verse mass of the (QQ) state (for the 3 S[ 8 ^ channel only): 



xap = \ 2 f Q) (B.2) 

where px (pj-) denotes the transverse momentum of the final-state gluon without 
(with) AP evolution and ^(qq) denotes the mass of the resonance. At high p?, 

p £P = XAp x pT (b.3) 

where pt is the transverse momentum of the resonance as generated by PYTHIA 
(i.e. without AP evolution), whereas for px < m (QQ) the effect becomes much 
less significant as it should be. Thus the interpolation between low and high 
Pt is smooth with the right asymptotic limits at both regimes. The above way 
to implement AP evolution may appear somewhat simple but it remains in the 
spirit of our whole analysis, i.e. using PYTHIA machinery whenever possible. 
In fact, it provides an energy depletion of the fragmenting gluon in agreement 
with previous work on charmonium hadropro duct ion |2(], In Figure 8 the 
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xap factor is plotted as a function of the transverse momentum of the resonance 
for the Tevatron event generation. 




1 9 2 C 
, (GeV) 



Figure 8: xap factor as a function of pt for Tevatron energies obtained from our generation. 



Moreover, in order to assess the effect of AP evolution on the fit parameters we 
show in Table 8 two numerical values for the relevant colour-octet MEs obtained 
from a best x 2 fit to Tevatron data [26f using the CTEQ4L PDF: (i) the first 
row corresponds to a generation without AP evolution; (ii) the second one does 



lT(l%3 



take it into account. Notice the increase of < O 
case w.r.t. AP off, but to a lesser extent than for charmonium W% . 



Si) > 



Mot 



in the latter 



Table 8: Colour-octet matrix elements (in units of 10 -3 GeV 3 ) from the best fit to CDF data at 
the Tevatron on prompt T(IS') production. The CTEQ4L PDF was used with AP evolution off 
and on respectively. 



ME: 


< Q? ,U "(W > L 


AP off 


70±15 


AP on 


77±13 



It is worth noting that the effect of the AP evolution on the shape of the 
differential cross section over the [1,20] GeV p^-range, though sizeable, is con- 
siderably less pronounced for bottomonium than for charmonium [12| because 
of the larger mass of the former. Nevertheless we can appreciate in Figure 9 
that the plot corresponding to AP evolution is noticeably steeper at moderate 
and high as could be expected. Let us finally remark that, although we can 
switch on/off AP evolution and initial-state radiation at will in the event gener- 



2(3 




Figure 9: Theoretical curves obtained from a fit using PYTHIA including the colour-octet mech- 
anism for prompt T(IS') production against CDF data at the Tevatron a) without AP evolution 
of the fragmenting gluon, b) with AP evolution of the fragmenting gluon. The CTEQ4L parton 
distribution function and rrn, = 4.88 GeV were employed in the fits; dotted line: CSM, dot-dashed 
line: 3 s[ 8 ^ contribution, solid line: all contributions. 



ation, both next-to-leading order effects have to be incorporated for a realistic 
description of the hadronic dynamics of the process. 
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Rapidity cut and azimuthal correlations 

In this appendix we show that the (systematic) uncertainty associated to the 
determination of the Feynman x of the interacting partons in our proposed 
method, is given by the upper rapidity cut yo on the resonance imposed to 
events, according to the expression: 

— = yo (Ci) 

X 

Indeed, assuming a gluon gluon scattering process into two final-state gluons, 
it is easy to see that any extra (longitudinal) rapidity amount Ay of any final- 
state parton, should be assigned to anyone of the two colliding partons, as a 
consequence of conservation of energy-momentum. (In this case the partonic 
reference frame would not longer coincide with the Lab frame.) 

On the other hand, a parton carrying a fraction x of the total hadron momentum 
has a (longitudinal) rapidity 

V = yhadron ~ log - (C.2) 
X 

where yhadron is the rapidity of the hadron in the Lab system. 

Differentiating both sides of Eq. (C.2) and setting Ay = y , one gets easily the 
expression (C.l). 

Let us observe that the rapidity cut \y\ < y binds us to a region of "allowed" 
transverse momentum, increasing with p T since x 2 s ~ 4p§,, and hence 

A P T (n i\ 
= yo (C.3) 

Pt 

This means that as the transverse momentum grows, the range compatible 
with the relative error, predetermined by choosing the value of yo, grows too. If 
this value is set very low, the precision on the Feynman x increases but the price 
to be paid is probably reducing too much the statistics. Conversely, allowing y 
to be too large, leads to larger statistics but spoiling the knowledge of x because 
of the uncertainty given by (C.l). As a compromise, we chose y = 0.25 which, 
however, could be varied depending on the size of the sample of collected events. 

In order to get an idea of the expected impact of the intrinsic kx on the topology 
of events, we show in Figure 10 several plots of the transverse momenta of the 
T(IS') resonance versus the recoiling jet. In the absence of any higher order 
QCD effect, events squeeze along the diagonal. However kr smearing spreads 
events over a larger area in the plot, spoiling somehow a naive picture of a 
back-to-back topology coming from a collinear approximation to leading order; 
Figure 10. a) corresponds to initial-state radiation activated in the PYTHIA 
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Figure 10: Plots of the jet transverse momentum versus the T(IS') resonance transverse momen- 
tum at LHC energy (parton/particle level simulation) using, from left to right: a) the PYTHIA 
algorithm to simulate initial-state radiation; using a gaussian smearing function with b) < >= 2 
GeV and c) < kx >= 3 GeV. The two straight lines indicate the allowed region according to the 
Pt uncertainty obtained from Eq. C.3 for yo = 0.25. 



generation following the model developed in |ZT]. Alternatively, Figures 10. b) 
and 10. c) show the effect of a gaussian spread of < kx >= 2 GeV and < kr >= 3 
GeV, respectively. The region inside the two straight lines corresponds to an 
uncertainty on px given by Eq. (C.3) for a rapidity value y = 0.25. Although 
at small and moderate pr (say, pr < 10 GeV) all plots essentially agree, at 
higher px the former one, corresponding to a full simulation of gluon emission 
in the initial-state performed by PYTHIA, displays much more events outside 
the accepted region. 



Table 9: Fraction (in %) of events inside the region defined by the two straight lines for different 
Pt lower cuts (in GeV) applied to the resonance, corresponding to Fig. 10. a), i.e. initial-radiation 
generated by PYTHIA. 



Pt cut-off: 


10 


20 


30 


40 


50 


% "inside" 


39 ±1 


38 ±3 


35 ±5 


38 ±9 


38 ± 15 



In Table 9 we show the fractions of events inside the allowed region between the 
two straight lines in the plot 10. a) (initial-state radiation on). We observe that 
about 40% of all events are "accepted", remaining practically constant above 
Pt = 10 GeV. Finally we conclude that such reduction factor (of the order 
of 40%) does not represent in itself a dramatic loss of statistics regarding our 
proposed method to probe the gluon density in protons. On the other hand, for 
the gaussian smearing, the situation is even much more optimistic. 

In Figures 11 we show the azimuthal A0 angle between the muon pair direction 
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Figure 11: Azimuthal angle between the recoiling jet direction (denned by the parent gluon momen- 
tum) and the dimuon direction from Y(IS') decays in the transverse plane, from left to right: a): 
Initial-state radiation activated in PYTHIA; b) Using gaussian smearing with < kr >= 2 GeV; c) 
The same with < kr >= 3 GeV. All plotted events were selected with T(15) transverse momentum 
greater than 10 GeV. 

(defining the direction of the fragmenting gluon into bottomonium) and the 
recoiling jet generated by the final-state gluon, for different values of the effec- 
tive kr, in correspondence with Figures 10. In Fig. 11. a) we used the PYTHIA 
algorithm for initial-state radiation, whereas in Figures 11. b) and 11. c) we used 
a smearing gaussian with < kr >= 2 GeV and < kr >= 3 GeV, respectively. 
As expected, again we realize the sizeable effect of the effective hp effect on 
the distribution, especially in the former case. Nevertheless, most events should 
display a clear enough back-to-back signature as regards the A0 variable (in 
addition to the p T balance), as indicated by the peak at 180 degrees in all plots 
of Figure 11. 
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